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emission of the photoelectrons and secondary electrons, 
the differential charging of the surface of the probe and 
a spatial potential distribution in the vicinity of the 
space probe. These effects and their importance are dis- 
cussed with consideration of the following plasma conditions: 
1) geostationary satellite orbits, 2) in the solar wind 
(HELIOS mission), and 3) in the ionosphere at an altitude of 
250 km (the projected OSV on Spacelab). A review of the 
fundamental models is given based on the kinetic theory 
of plasma (Vlasov-Maxwell system). 
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INTERACTIONS BETWEEN SATELLITES AND PLASMA 


U. Isensee, W. Lehr, H. Maassberg 
Applied. Geophysics, Institute for Applied Physics, 

Technlsche Hochschule, Darmstadt 

1. Introduction /305* 

A real problem arises In the measurement of plasma from 
satellites: the surrounding plasma is disturbed by the satellite 

Itself and Its surface is electrically charged. For the evalua- 
tion of the plasma measurements, the satellite's potential and the 
disturbance must be known in order to be able to determine the 
undisturbed plasma. Corresponding to the schematic presentation 
in Fig. 1, space charges around the satellites are produced by 
the following effects: 

a) The relative streaming plasma is faded-out by the satel- 
lites. Due to their very high thermal speed, the electrons 
refill the wake very strongly, so that an extensive nega- 
tive space charge is formed. 

b) Photoelectrons, which can form a cloud of negative space 
charges in front of the satellite, are emitted from the 
parts of the probe's surface which are irradiated by the 7306 
Sun. 

c) Secondary electrons can be emitted by the reflection of 
high energy plasma particles. 

d) Plasma electrons are backscattered on the probe's surface. 

The spectrum of the backscattered electrons contains both 
elastic and inelastic stray components. 

These interactions cannot be experimentally simulated in 
plasma chambers. Even the chamber's walls would emit secondary 


^Numbers in the margin indicate pagination in the foreign text. 
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electrons and backscatter 
the plasma electrons — to 
a much greater degree than 
the small simulation probe 
itself. Furthermore, the 
necessary thermal isotropy 
and homogeneity of the 
plasma's electron components 
cannot be attained. 


Fig. 1. A schematic representation 
of the space charges around a 
satellite . 


Key: 1. Streaming plasma; 2. Sat- 
ellite; 3. Sunlight; 4, Photoelec- 
trons; 5* Wake; 6. Ions; 7* Elec- 
trons . 


The surface is electric- 
ally charged by the plasma 
stream on the probe and by 
the streams of photo- and 
secondary electrons moving 
away from the probe. The 

floating potential is obtained from the flux balance (the disappear- 
ance of the joint stream). The knowledge of this floating potential 
is, for example, necessary for the interpretation of the measured 
electron spectrum in the HELIOS mission. The properties of the 
surface material are included in the flux balance, which for dif- 
ferent parts of the surface results in the general differential 
potential (differential charging) . Potential differences of around 
10 kV between the satellite's various parts can be obtained for 
geostationary satellites, which has already led to electronic fail- 
ures in the satellite. The SCATHA (Spacecraft Charging at High 
Altitude) program of the USAF should experimentally study these 
charging effects . 


Space plasma is very homogeneous in comparison to that in 
a plasma chamber, so that it is especially suitable for the study 
of wave phenomena. The potential distribution around the probe, 
which is induced by the undisturbed plasma, can stimulate the 
waves themselves (instabilities). There is, however, only little 
information available about this type of instability taking into 
consideration the magnetic field. 
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The charging 0/ a probe by particle fluxes 


Pig. 2. A schematic representation of 
the flux balance on an element of the 
probe’s surface for plasma conditions in 
a solar wind. The scattering of the 
photoelectrons on the other parts of 
the satellite’s surface is essential. 

The potential minimum in front of the 
probe represented here is obtained 
only in a solar wind, in other cases 
different surface potentials cause a 
similar effect . 


The need, therefore, 
arises to investigate 
the Interaction of 
satellites with space 
plasma in order to reach 
the first level of a 
universal probe theory 
valid under conditions 
of outer space. If 
the plasma in a strong- 
ly disturbed region is 
not in thermal balance 
(great free path), the 
necessary theoretical 
treatment can ensue 
only on the basis of 
the statistical kinetic 
theory . 


2. Basic Equations 

2.1. The Vlasov-Poisson System 

The following approximations are accomplished very well for 
all of the plasma conditions discussed in this article: 


a) The interactions between satellites and plasma should be 
stationary. The possible stimulation of instabilities 
can be ignored. 

b) The plasma can be treated in a very good approximation 
as smooth (Vlasov plasma). 

c) The influence of the magnetic field can be ignored because 
of the very large gyro radii relative to the satellite's 
characteristic dimensions. This supposition is only 
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condlt 3 onally valid for the ionospheric conditions. 


Consequently, the Interaction of the satellite with the surrounding 
plasma is to be described by the Vlasov-Pod sson system: 


The Vlasov formula: 

for p = i, e the plasma particle types 


The Poisson formula: 

OO 

a<k*) - d5v 

r -do r 


( 1 ) 


( 2 ) 


The flux balance formula: 


i(*-X s ) - JZ % [If Z f u ( 2.Z) d 5 v 


hK 1 ' 

with x on the probe's surface; 




(3) 


The flux's normal component disappears for a surface element x_ 

— s 

with Insulated material. 

Z \\ A(x s )-*l - 0 

“ ? 

The total flux on all the conductive surface parts F^ (which are 

s 

electrically connected) disappears. An electrical model of the 
satellite is necessary here (with a time-dependent spin). 


The formulas contain the following symbols: 
x, v location and velocity coordinates 

q , m charge and mass of the p-type particles (p = i, e) 

H r* 

f (x, v) distribution function in the 6-dimensional phase space 

H 

$(x) electrical potential 

j_(x) electrical current density 

The following limit conditions are valid for a definite 
solution of the Vlasov-Poisson system: 
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f i (x, v) = 0 with v*n >_ 0 


All the Ions should be neutralized on the surface. 


f e ( 2S B *£> 




* * 




With, v-n^ 0 


Photo Secondary backscattered electrons 
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in general for electrons isotropic Maxwell distribution, for ions, 
Maxwell distribution shifted by the relative velocity; 

lim - 0 

|Kt-*oo 


The other boundary condition for the potential <j> on the probe 

<Kx e ) = <f> (x) must be determined from the flux balance equation (3). 
— fcj b — s 


This Vlasov-Poisson system with boundary conditions can be 
solved with the characteristic methods: the distribution function 
f (x, v) is constant in a particle's orbit (a characteristic of /308 

r* 

the Vlasov formula). The distribution function can be determined 
by orbital tracking up to the limit value. The combined particles 
(trapped particle orbits) can be ignored here. 

2,2. Iteration of the Vlasov-Polsson System 

A general solution technique of the integration of the Vlasov- 
Poisson system will be presented below in order to clarify the 
causal relationship of the system and of the limit value. An 
appropriate selection of the initial potential <j>° (x) is presumed, 
the upper index indicates the corresponding stage of the iteration. 

Potential ^(x) 

1. The calculation of the distribution function f^(x, v) in 
the potential <|> (x) is carried out by the particle orbit tracking 
up to the limit value on the probe's surface or in the undisturbed 
plasma. 
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(In general, the distribution function for the electrons 
emerging on the probe's surface is dependent on the distribution 
of the arising particles (secondary electrons and backscattering) , 
so that the distribution function of the emerging electrons must 
first be calculated by a further Iteration.) 

2a. The particle density n^(x) is calculated by integration 
over all the possible orbits: 

n k (x) ° 

2b. The streams on the satellite's surface are calculated: 
Prom the flux balance (3), then follows the boundary condition 
for the potential on the surface x g : <J> k+1 (x. s ). The parts of the 
probe's surface are the capacities against each other and against 
the undisturbed plasma). 

3. The solution of the Poisson formula (2): from this 

results the potential or the next iteration step 

Potential <t> k+1 (i 

In general this iteration procedure is convergent. The self- 
consistent potential and density distributions obtained thus are 
then the solutions of the Vlasov-Poisson system (1-3). Calcula- 
tions of this type have been performed numerically with great 
simplifications (for example, by I. Katz et al . , 1977, and L.W. 
Parker, 1977)- 

2.3. Numerical Plasma Simulation 

A completely different solution procedure is numerical plasma 
simulation, in which the motion equation of very many particles 
is simultaneously solved (for example, by Isensee, 1977). The 
plasma particles are represented as discrete particles of very 
great mass and charge because of the numerical capacity restric- 
tion according to the particle-in-cell method (Morse, 1970). 

This is performed on the outer limit of the simulation region 


and on the probe's surface which corresponds to the boundary value 
of the distribution function. This simulation particle moves for 
a long period of time Ir. the potential, which is then calculated 
anew through the tabulation of the particle densities and the 
solution of Poisson equation (2). Each simulation particle is 
observed until Its orbit either leaves the simulation region or 
reaches the probe's surface and contributes to the flux balance - 
The floating potential is thus regulated so that the common stream 
on an element of the surface disappears. The simulation Is per- 
formed until the potential no longer changes substantially. Sta- 
tistical fluctuations In the particle densities and the potentials 
result from the small number of simulation particles (around 10^). 
Another numerical problem Is the very small spatial resolution near 
the probe's surface (see also D.W. Parker, 1977). In this form 
the numerical plasma simulation is only possible if, on the one 
hand, the simulation region is several Debye lengths long, and, 
on the other hand, a sufficient number of simulation particles are 
inside a Debye sphere — this is accomplished directly with the 
plasma conditions of solar winds (see Section 3.2). The stationary 
potential and density distributions are then the self-consistent 
solutions of the Vlasov-Poisson system. 

3 . Models Under Different Plasma Conditions 


3.1. Plasma Conditions 

The table on the following page cites the typical plasma 
parameters for the conditions in solar winds, for geostationary 
orbit and in the ionosphere, respectively. 

A highly thermal plasma with very slight density is in a 
geostationary orbit; in the ionosphere it is cold, but very dense. 
In solar winds somewhat more moderate conditions are encountered. 
Apart from the electrical components the effect of the magnetic 
field is negligible in the ionosphere. 
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TABLE X. THE TYPICAL PLASMA CONDITIONS FOR A SOLAR WIND, GEOSTA- 
TIONARY ORBIT AND MIDDLE IONOSPHERE (PER DAY). 


Plasma parameter : 


Theri al energy 

e 

kT. 

Plasma density n Q : 

Relative velocity V R 
of plasma to probe: 

Power numbers S e,i : s p" v R/ v thti 

Dpbye length ^ 

(L 3 £fc t : Probe dimensions) ‘ D 

Effect of magnetic field 
(r^,: Gyro radius) 

Photoelectrons kT Ph 

n Ph 


Solar wi n( j 

O.J - i.o 

Geostationary 
Orbit ( 6.6 R e ) 

1 Ionosphere 
1 250 kc 

10 - 50 eV 

5-50 keV 

I 0,2 eV' 

10 - 20 eV 

5 - 50 keV 

1 0.1 eV 

10 - 500 cm " 5 

0.2 - 5 cra - ^ 

| 2 - 10 5 cm " 5 

500 km/eec 

50 km/sec 

7 - B km/sec 

s e « 1 « Si 

s e « S ± « 1 

s e « i« s i 

~ L . 
D aot 

>> L s 8 t 

*D « L S et 

r r » L . 
« sot 

r G» L B 6 t 

r Gi>> L S et 
r Ge <<: L eet 

1 - 5 ov 

1 - 5 eV 

1 - 5 eV 

0 ? - 10 4 cm“5 1 

JLO* - 10 4 cm “ 5 

10 ^ - 10 4 cm~^ 


3.2. Solar Winds 

With the aid of numerical plasma simulations the potential 
distributions and the densities around a two-dimensional model 
probe (corresponding to the HELIOS satellites) were calculated 
by Isensee (1977 and 1978). In Fig. 4, the extensive negative 
potential structure can be seen through the fade-out of the solar 
winds’ protons in the wake of the probe. A negative potential 
barrier, which is caused by the very dense photoelectric clouds, 
is found in front of the probe. All the electrons are presented 
in Fig. 3, while only the photoelectrons are in Fig. 5- The 
solar wind electrons, although faded-out by the probe, are only 
a little disturbed by the potential. They, like the protons, 
form a constant background. Secondary electron emission and 
backscattering are negligible under these plasma conditions. 

In front of the illuminated parts of the surface (the Sun shines 
from the left), the photoelectron density is very great and the 
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QUALITY 

inserted line corre- 
sponds to the position 
of the half-disturbed 
plasma density. The 
photoelectrons are 
scattered on the poten- 
tial in the wake and 
cannoc penetrate it. 

The solar winds' 


Pig. 3. The simulation particle in the 
numerical plasma simulation for a rec- 
tangular probe model with plasma condi- 
tions corresponding to those of a solar 
wind. Only the electrons are represen- 
ted here. Left in front of the probe 
is the dense photoelectron clo'..; cor- 
responding to Pigs. 4 and b. The Sun 
shines from the left (2-dimensional 
simulation model according to Isensee, 
1978). The distance from the Sun: 

R « 0.3 AU. 


protons are only 
negligibly disturbed 
by the probe. Because 
| e$ | << kT 4 << Ej in 
(Ej ln •* proton slip 
energy;, the ion den- 
sity behind the probe 
can be described by 
the neutron approxima- 


tion. On the other 

If i n 

hand, £* << |e<f>| << kT p is valid; a local Boltzmann distribution 

with geometric fade-out by the probe can therefore be taken as a 


good approximation of the solution of the Vlasov equation (1) for 
the 3olar wind electrons. Therefore the Vlasov equation (1) must 


still be solved only for the exiting photoelectrons: this is 
certainly not possible due to the strong nonlinear aspect of the 
Vlasov-Poisson system for realistic conditions. 


For the wake, where the photoelectron density does not nlay 
an important role (Debye lengths in the probe's dimensions), the 
potential cah, however, be analytically calculated (Fig. 6). 

For shorter distances from the Sun, the wake's structure is sig- 
nificantly more pronounced. The floating potential in Fig. 6 
stems from a spherically symmetrical (analytical) photoelectron 
model (Maassberg, 1978b). Here, because of the symmetry of the 
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potential in tht vicinity ot a probe 


n.78cm-l r,.0 25 t0 e K w.. 500 km/*«c 


Fig. The potential in the vicinity of a two- 
dimensional model probe in a perspective repre- 
sentation (w * slip velocity of the solar wind, 
plasma parameter at around 0.3 AU. According 
to Isensee, 1978). Analogous to Firs. 3 and 5. 



Fig. 5 . Numerical plasma simulation. Photoelec- 
trons in the vicinity of a probe in a solar wind 
(analogous to Figs. 3 and 4). 


Vlasov-Poisson system, ( 1 ) and ( 2 ) and the flux balance equation 
can be approximated. 


Under the conditions of a solar wind, the strong nonlinear 
character of the Vlasov-Poisson system is the main difficulty for 
describing the interactions of a probe with the surrounding 
plasma. Model calculations of this kind are necessary for 
interpreting the electron spectra, which are partially strongly 
disturbed by the potential near the probe, measured by the HELIOS 
mission. 
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3.3. Geostation- 
ary Orbit 

The plasma in 
this region (radiation 
belt) has a much higher 
temperature than in a 
solar wind, but a lower 
density. The essen- 
tially greater electron 
flux negatively 
charges the probe's 
surface quite strongly, 
so that no photoelec- 
tron cloud can be 
formed in front of 
the probe. Since 
according to Table 1 
the Debye length is 
very large In compari- 
son with the charac- 
teristic dimensions 
of the probe, all of 
the space charges vis-a-vis the surface be ignored. The float- 
ing potential is practically unshielded. Consequently, instead of 
the Foisson equation (2) only the Laplace equation must be solved 
for tne potential: A<|>(x) * 0. In order to calculate the floating 
potential as the boundary value, the flux balance equation (3) must 
be integrated self-consistently with both the Vlasov equations (1). 
Photoelectrons, secondary electrons and backscattering near the 
plasma stream on the surface determine the potential structure 
around the probe. These streams are self-consistently calculated 
by Katz et al. (1977) for different parts of the surface by 
orbital tracking according to the iteration process in Section 2.2. 
Certainly the secondary electron emission and backscattering are 
not taken into consideration. According to the model calculations 



Fig. 6. A perspective representation of the 
potential in the wake of a sphere (radius 
1.2 m) with plasma conditions for R * 0.3 AU 
(above) and R = 1 AU (below) according to 
Maassberg, 1978a. 
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of Prokopenko and Lafram- 
bolse (1977) , however, 
these effects must be 
considered . 


The potential dlstrl- /312 
butlon around a square 
block with different sur- 
face materials Is pre- 
sented in Fig. 7. The 
cross-hatched region 
represents metal, the rest 
should be covered with a 
dielectric film (solar 
cells). Since the potential 
like that of a monopole in 
the center of the model probe is related to the distance from the 
surface, the effect of the different surface potentials subsides. 

Potential differences of several kV appear, however, on the surface, 
the pbotoelectrons cause e somewhat higher potential on the illumina- 
ted, insulated parts of the surface. The chief problem under these 
conditions is the self-consistent solution of the Vlasov-Laplace 
system in connection with the flux balance equation for the indi- 
vidual parts of tne satellite's surface with different materials. 

The properties of the materials with regard to the work function 
and photo- and secondary electron spectra have been recently 
investigated (for example, Prokqpenko and Laframboise, 1977). 

3_j£> The Ionosphere 

The iorosphere is represented fully contrary to the plasma 
conditions for a geostationary orbit: high plasma density and 
low temperatures. The Debye lengrh is very short in comparison 
with the probe's characteristic dimensions, so that the floating 
potential is very quickly shielded. The photo- and secondary 
electron emissions as well as the backscattering are negligible 
because of the high density and low temperatures. The thermal 



Fig, 7. Equipotential lines around a 
model satellite with conductive (cross- 
hatched) and nonconduct ive surfaces. 

The Sun shines from the right. (Acc- 
ording to Katz et al., 1977.) 



velocity of the plasma electrons Is very great In comparison with 
the slip velocity of the ions; the surface is hence electrically 
charged quite negatively, so that only a very small part of the 
electrons reach the probe. Furthermore, the electrons' gyro-radius 
is very small relative to the probe's characteristic dimensions 
(adiabatic motion of the electrons along the line of the magnetic 
field in the region of the probe). Therefore, the electron fade- 
out by the probe relative to the density is negligible. The 
potentials around the probe are negative everywhere, therefore the 
electron density can be approximated very well by a Boltzmann 
factor: n e (x) = n 0 exp(e<J>(x)/kT e ) . The Vlasov equation (1) for the 
electrons need not be solved any longer. For smaller probes 
( <1 . 0 m) the Influence of the magnetic field on the ions' motion 
can also be Ignored, with the plasma experiment planned within 
the framework of the OSV (Orbital Flux Experimental Station on 
the Spacelab) the probe to oe measured is in this range. With the 
exception of the region directly behind the probe (a much closer 
wake), the condition of quasineutrality can be taken because of 
the very small Debye length in relation to the probe's dimensions 
and therefore also in relation to the characteristic changes of 
the potential in the wake: n g (x) = n^(x). Therefore in order to 
describe the wake structure under conditions of the ionosphere 
only the Vlasove equation (1) can be solved for the ions. The 
electron density is approximated by a Boltzmann factor and the 
Poisson equation (2) is solved by the conditions of quasineutrality. 

The floating potential has no influence on the wake structure. 

3.5. The Influence of Surface Materials /312 

In the previously described models, the interaction of the 
plasma with a satellite relative to the potential distribution in 
the plasma has stood in the foreground. In the following section 
a simulation model, following Soop (1972), is presented in which 
the external plasma is ignored (for example in the solar wind at 
the far outer half of the Earth's orbit). Photoelectrons, which 
return to the probe because of the surface's positive potential 
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(the energy spectrum of the photoelectrons must be restricted to 
the top), are emitted from the probe's surface. The photoelectron 
clouds in front of the surface are shown in Pig. 8. The parts of 
the surface lying in the shadow (the picture on the right) have a 
small potential also on the front side in the sunlight, so that 
the photoelectron clouds are concentrated more in front of the 
illuminated parts. In the case of the conductive surfaces (left), 
the total probe is surrounded by photoelectrons, v/hose density in 
front of the illuminated parts of the surface is only somewhat 
larger . 



Pig. 8. The photoelectron clouds in front of a 
sphere ignoring the surrounding plasma. In the 
left picture is a conductive surface., in the right, 
an insulating one is simulated. The Sun shines 
from the right. Numerical plasma simulation 
according to Soop (1972). 


The surfaces of real satellites consist of both conductive and 
non-conductive materia] s (solar cells). The charging effect and 
the potential structure around the satellite are, therefore, very 
complex and only barely comprehensible with the models used 
heretofore. The future direction of work in the field of probe 
theory in outer space conditions will depend on the information 
obtained from satellite measurements (plasma spectra) being 
understood on the basis of the models discussed here and the 
development of more realistic models. Similarly the results to be 
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expected from the SCATHA program of the USAF are of special sig- 
nificance for the direction of modeling in the field of geosta- 
tionary orbits and the plasma experiments planned within the 
framework of OSV on the Spacelab for the calculation of the 
wake structure under conditions of the ionosphere. 
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